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Abstract 

Static lookup tables as often used in the position calculation electronics 
in position sensitive detectors suffer from the well known problem that 
the propagation of digitization errors in the division leads to unequal 
efficiencies for the different output channels. In this paper we present a 
solution to the construction of such a lookup table which is optimal in 
the sense that the least possible efficiency errors are committed in the 
class of monotone lookup tables. The obtained resolution is comparable 
to what is obtained using the standard implementation of the fixed point 
division. After a theoretical explanation, we present some experimental 
results confirming our claim. 



1 Introduction 

Resistive charge division is a well-known technique for determining the position 
of the impact of an ionizing particle along a resistive electrode read out by 
charge amplifiers on both sides. If X is the charge measured at side A and Y 
is the charge at side B, then the position can ideally be calculated from the 
dimensionless number: 

P=^— r (1) 
X + Y V ' 
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Several electronic techniques exist to calculate the relative position P from the 
electronic signals X and Y. Although analog treatment has been used in the past 
with success, the progress in ADC converters and digital circuitry promotes an 
all-digital solution. Moreover, the existence of large EPROM memories invites 
the use of an extremely simple design, which is often used in correspondence 
with position sensitive detectors: a lookup table which associates which each 
couple (X, Y) the position P. In practice, we've build a circuit with two 11-bit 
inputs for X and Y and one 8-bit output for P. 



2 Digitization Problem 

Although the task of filling in the lookup table may seem simple (use the dig- 
itized values of X and Y, and calculate P, rounded to 8-bit accuracy) it leads 
to an unequal efficiency of each of the 256 channels due to digitization effects. 
Perfectly uniform irradiation would then lead to unequal intensities observed in 
the different channels. In the rest of this paper, we will refer to this assignment 
as the 'standard position calculation'. The quantization effect and the reasons 
why they are sometimes a problem are well-described in pQ. In this paper we set 
out to find the best possible lookup table which reduces this effect maximally, 
without compromising the resolution. 

The effect can be explained as follows: a true value X will be represented by 
its digitized counterpart Xd and a true value Y will be represented by Yd- So 
a digital couple (Xd, Yd) represents a square in the (X, Y) plane with center of 
gravity (Xd + 1/2, Yd + 1/2). The standard way of assigning this point (and so 
all couples in the square we're considering) to a position bin is by calculating: 

p , = X d + l/2 

Xd + Yd + 1 V ' 

and after digitization: 

Here, N out is the number of bits in the output (for 256 position channels, it is 
equal to 8). The problem arises because this assignment maps the whole square 
of (X, Y) values to the value Pd as determined by the exact position of its center, 
and the distribution of these center points is not uniform over the 256 position 
channels. The probability for a point (X, Y) to be solicited depends on (X + Y) 
through the pulse height spectrum of the incident radiation. 

We make one fundamental hypothesis: we assume that the pulse height 
spectrum of the captured radiation is independent of position. So we're actually 
considering an event distribution which factorizes in the variables E and P: 

da = spcc(E)ima,ge(P)dPdE (4) 



We know spec(-E) (the pulse height spectrum, assumed independent of the po- 
sition), and we want to find out image(P) from a(E,P). However, technically, 
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we're working in the variables X and Y, so we have: 

da = p(X, Y)dXdY (5) 

So our task is to extract image(P) from p(X,Y). Working directly in the dig- 
itized variables (but analytically continuing them...), we have the following co- 
ordinate transformation: 

E = X+Y+l 

P = X + 1/2 (6) 
X+Y+l K ' 

The Jacobean of the transformation d(P, E) /d(X, Y) = 1/(X + Y + 1) = 1/E, 
so we obtain: 

p(xy) = !P^) image(p) (?) 

This suggests that points which are geometrically uniformly distributed in the 
X, Y plane (as our digitized (Xd, Yd) couples are) have a weight equal to spec(E) / E 
in the image. It would hence be logical that each bin of the image (each digi- 
tized slice of the variable P) contains as much as possible an equal total weight 
of spec(E) / E over each of the points (Xd, Yd) assigned to the bin in the lookup 
table. Of course there is the constraint that the P'-value of each couple should 
still be close to the nominal P-value of the chosen bin (otherwise the resolution 
of our image will suffer seriously). We propose a method which satisfies both 
requirements. 



3 Proposed method 

There are already several successful techniques that are proposed to tackle this 
problem: 0| proposes to assign the (X,Y) point partly to different P-bins ; 3 
extends the effective word length over which the division is worked out. 0] 
partly uses a lookup table to invert (X+Y) and uses a processor to calculate 
the final position. However, all these techniques change the basic electronic 
architecture and do not use a simple lookup table. |S] and jS] have written about 
half-analog methods, using logarithmic amplifiers. |o| has moreover proposed a 
more classical method of correcting the image off line with efficiency correction 
coefficients which can be calculated from the pulse height distribution. He points 
out that this gives rise to increased statistical errors. yQ does give an improved 
way of calculating a lookup table, but the method is sub-optimal, although 
probably sufficient for many applications. 

We set out to find a mathematically optimal solution of constructing a lookup 
table. An optimal solution can eventually be used as a benchmark to compare 
to other techniques. 

We first create a list L of all possible pairs (Xd, Yd) of which there are, in 
our case, 2 22 = 4194304. We define a function xy, that maps the set of integers 
{1..2 22 } onto the set of couples {0..2 11 — l} 2 in lexical order: 



xy(i) = (X d ,Y d ) 



(8) 
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Next, ordering the couples according to their P' value given in equation[21comes 
down to the definition of a permutation p of 2 22 elements, defined as follows: 

P'{xy(p(i))) < P'(xy(p(j))) &i<j (9) 

The list L' of the pairs xy(p{i)) when i is running from 1 to 2 22 is the list of 
position-ordered couples. We assign a (non normalized) weight to each element 
in this list in the following way: 

Wl = spec(X d + Y d + l)/(X d + Y d + l) (10) 

with 

(X d ,Y d ) = xy(p(i)) (11) 

Taking the normalized, cumulative sum in the order of the list L', we finally 
obtain: 

n = pi^ (12) 

The new mapping is now defined as: 

xy(j) i-> [256r p -i (j) ] (13) 

Here we used the inverse permutation p~ l , to obtain a final mapping in lexical 
order, which is necessary to practically program the EPROM memory. In this 
way, each of the bins has very close to the same summed weight of points in 
it. We also respect the order of the value P' so that a couple with a higher 
P' value can never be assigned to a lower bin than a couple with a lower P' 
value. Within the family of all possible lookup tables which respect the order of 
P' , the proposed method gives us the optimal solution concerning the flatness of 
the efficiency, by construction. Indeed, the construction here proposed makes 
the efficiencies (the sum of the weights) equal, up to one single effect, that is: in 
the order of P', one has to decide whether a last point will still be attributed to 
the lower bin (and in that case, the extra weight of this single point might give 
us a total sum slightly above 1/2 ), or we might decide to attribute it to the 
next bin so that the lower bin has a total sum slightly under 1/2 N . Any lookup 
table construction respecting the order of P' will have to make this decision and 
hence at least this error, but in our case, it is the only error. Moreover, apart 
from this decision on these pivotal points, the proposed lookup table is unique. 
Hence any other lookup table will give worse efficiency errors. Although quite 
computing intensive, resolving for the permutation of 4 million pairs of numbers 
is now within reach of a modern personal computer and can be accomplished in 
a matter of minutes. 

Given that we have now the optimal solution concerning efficiency flatness, 
we will have to investigate what happened to the position resolution. It turns 
out that the error in position thus committed is not much larger than the normal 
quantization error of the standard method (which is optimal with respect to the 
position resolution). This is understandable, because in the limit of very high 
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binning in X and Y (for the same resolution in P), that is, in the continuum limit, 
the difference between the standard method and our newly proposed technique 
vanishes. Indeed, if we work in the real variables (X,Y), limited to the square 
[0, l] 2 , let us define digitized variables: 

(14) 

with N the very wide word length of the converters. The accumulated non 
normalized weight in output channel n (one of the 256), using our new definition 
of weight, but using the standard assignment, can be approximated very well 
by an integral: 

$>i=i / / spec{E )/E— — (15) 

ien Jo<E<E max Jn-l<256P<n e e 

E w -i/ spec(E)dE (16) 

We now see that the standard assignment already makes the weights accumu- 
lated in each position bin equal, so we would obtain exactly the same assignment 
using our new technique (because the lookup table doing this is unique), what- 
ever pulse height spectrum is used. It also means that if we can estimate the loss 
in resolution (which will turn out to be insignificant) in a low-bit word example, 
this will give us an upper limit on the resolution loss. We will study a Monte 
Carlo example to have an idea of what happens. 



X d = 



X 




Y 




;Y d = 




e 


e 



4 Monte Carlo simulation. 

In order to illustrate the problem and its proposed solution, we simulate the 
behavior in the case of a low-bit example: we suppose that the incoming signals 
X and Y are digitized on 6 bits, and that we calculate the position on 6 bits (64 
positions). We take as a spectrum an inverted Raleigh distribution, leading to 
the probability density function: 

spec(£) = 0.00781(63 - £) e -o.oo39(63-iJ) 2 

The precise choice of this distribution doesn't really matter, it just grossly 
looks like a true thermal neutron spectrum (using He-3 as a converter gas and 
in the proportional region) and is mathematically well defined. The proposed 
"spectrum" is displayed in figure ^ We simulated 10 6 points with a uniform 
position distribution and an energy distribution drawn from the above spectrum. 
Introducing the digitization according to a 6 bit scale for X and Y, we observe, 
using the standard position calculation, a severe digitization influence (with 
spikes of the order of 25 %) on the uniformity of the response of the system, as 
displayed in figure 2, full line. 
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Figure 1: Artificial distribution that will serve as 'pulse height spectrum' 



Using the newly proposed technique, we obtain the distribution of the cal- 
culated position in figure |2 dashed line. The uniformity is clearly superior in 
the this case. The price to pay is a very small decrease in resolution. Let us 
investigate how severe it is. The quantization error in the case of the classical 
calculation is of course close to a uniform distribution with a standard devi- 
ation of 1/V12 = 0.28867. Using the Monte Carlo data, we find a standard 
deviation of 0.292. The standard deviation of the error using the new mapping 
is measured to be 0.298 when using a uniform weight over all couples (X,Y). 
The two distributions of the position errors are shown in figure 01 We hence 
observe a relative loss in resolution of no more than about 2%. Essentially, the 
resolutions can be said to be equivalent, because this tiny decrease in resolution 
shouldn't affect the spatial resolution of the overall system, which shouldn't be 
determined by the number of output bins, but by physical processes and noise 
limiting the intrinsic resolution of the detector. 

If we quantify the non-uniformity by calculating the standard deviation of 
the bin contents (using our sample of 1000000 events) we find the residual 
standard deviation before quantization is 133.3 counts (per cell) where we expect 
125. Using the standard way of calculating positions, we have a non-uniformity 
of 1686.7 counts, and this is reduced to 388.5 counts using the newly proposed 
technique. 

However, in order for this to hold it is important that the true spectrum in 
the data matches closely the spectrum used in the design of the lookup table. 
If we shrink the spectrum of the data by 5% (applying a factor 0.95 to the orig- 
inal spectrum used to construct the table), the good uniformity of the position 
histogram is partly gone. The standard technique gives us a non-uniformity 
of 1902.3 counts while the new technique obtains 794.2 counts. A shrinking of 
10% leads to a standard non-uniformity of 2236.3, while the proposed technique 
obtains 1424.0 counts. Although still better, this indicates the need to match 
closely the spectrum used in the composition of the table to the actual spec- 
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Figure 2: Result of the standard and new assignment of events to position 
channels. 




Figure 3: Distribution of position errors for the standard technique and the new 
technique. 



Optimal Lookup Table in Charge Division 



8 



normalized 



intensity 



0.002 




2000V 



1950V 



0.0015 



0.0005 



0.001 




channel 



500 



1000 



1500 



2000 



number 



Figure 4: Normalized Pulse height spectrum of the detector at 2000V and at 



trum that will be used in order to take the full advantage of the method. For 
proportional counters used in thermal neutron detection, this is a reasonable 
requirement as the spectrum is very stable (and also rather broad). On the 
other hand, given the very small word lengths in this example, this simulation 
is a particularly severe test for our method and in practice the sensitivity to a 
change in spectrum is smaller, as will be shown in the next part. 

5 Experimental verification 

Using an Am-Be thermal neutron source, a 40 cm long position-sensitive He-3 
neutron detector available from Reuter-Stokes of the type RS-P4-0814-2, and 
two 3V/pC amplifiers with an overall gaussian shaping time of 1.4 microseconds, 
we bias the detector at 2000V. This uses almost the full dynamics available (2048 
channels for 5 V) without saturation as illustrated in figure 0] 

In order to have very high statistics in a reasonable amount of time, we put 
the source very close to the detector; we will have no uniform illumination, but 
the smooth bell form image will indicate local fluctuations in counting efficiency 
also clearly. Using the standard algorithm and the new technique, we obtain 
different images as displayed in figure[S]with the dash-dotted and the full curve, 
respectively. 

The digitization glitches are clearly visible in the first one and buried in 
the statistical noise in the second. (Note that the position of the source is 
slightly different, as it was a manual positioning in between experiments). Let 
us now look at the sensitivity of the spectrum: we change the voltage from 
2000 V to 1950 V and use the same lookup table (constructed for use with 
the 2000 V spectrum). Notice that the applied spectrum (see figure 0] dotted 
curve) now has changed by almost 30%. When looking at the resulting image in 



1950 V. 
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Figure 5: Images obtained with the standard method and with the new method 
using the spectrum at 2000 V. The full curve shows the new image at 2000V, 
and the dashed curve shows the new image at 1950 V. The dash-dotted curve 
gives us the image using the standard method (the number of counts has been 
divided by 20 for this last image in order to fit on the figure). Because of the 
different manual positioning of the neutron source, the standard and new images 
do not coincide perfectly. 

figure [3] (dotted curve), we notice that a small digitization error appears around 
channel 128, but it is still much smaller than if we would have used the standard 
algorithm. 

We have hence experimentally illustrated that the glitches in efficiency dis- 
appear when applying our new technique. We've also demonstrated a certain 
robustness against differences between the pulse height spectrum used to con- 
struct the lookup table, and the pulse height spectrum of the detector using the 
lookup table. 
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